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ABSTRACT 

We use high-resolution hydrodynamical simulations to investigate the spatial correlation be- 
tween weak (TVhi < 10^^ cm~^) Lya absorbers and gas-rich galaxies in the local uni- 
verse. We confirm that Lya absorbers are preferentially expected near gas-rich galaxies 
and that the degree of correlation increases with the column density of the absorber The 
real-space galaxy auto-correlation is stronger than the cross-correlation (correlation lengths 
''"o.gg — 3.1 ± 0.1 Mpc h^^ and rn.nn — 1.4 ± 0.1 Mpc h^^, respectively), in contrast with 
the recent results of (lRvan-Webej2006l RW06), and the auto-correlation of absorbers is very 
weak. These results are robust to the presence of strong galactic winds in the hydrodynamical 
simulations. In redshift-space a further mismatch arises since at small separations the distor- 
tion pattern of the simulated galaxy-absorber cross-correlation function is different from the 
one measured by RW06. However, when sampling the intergalactic medium along a limited 
number of lines-of-sight, as in the real data, uncertainties in the cross correlation estimates 
are large enough to account for these discrepancies. Our analysis suggests that the statistical 
significance of difference between the cross-correlation and auto-correlation signal in current 
datasets is ^ 1-cr only. 

Key words: intergalactic medium, quasars: absorption lines, galaxies: statistics, large-scale 
structure of universe 



1 INTRODUCTION 

Understanding the interplay between galaxies and the intergalac- 
tic medium (IGM) is a fundamental cosmological problem. On one 
side the IGM acts as reservoir of gas that cools down in the poten- 
tial wells of dark matter haloes and forms galaxies and stars. On the 
other side the IGM is a sink that records, over a large fraction of the 
cosmic time, the crucial thermal and chemo-dynamical processes 
related to galaxy formation. Significant progress has been made in 
the last few years thanks to high-resolution spectroscopic data from 
quasar (QSO) lines-of-sight and imaging of QSO fields that has 
been performed by several groups. The properties of Lya and metal 
absorption lines in the high redshift universe have been cross- 
correlated w it h thos e of the galaxies (e.g. Adelberger et al. (200j); 
Nestor et al.l j2007h; lBouch6 et al.l ( 12007^ : .Schaye et al.. C20d7h : 
Churchill et al. to shed light on the physical state of the 

IGM around them and possibly on the still poorly understood feed- 
back mechanisms. Among all the possible elements in various ion- 
ization stages hydrogen is the most abundant and thus has been 
widely studied by the scientific community. The analysis of the sta- 
tistical properties of Lya lines and of the transmitted flux shows 
that the neutral hydrogen in the high-redshift universe is embed- 



ded in the filamentary cosmic web that traces faithfully, at least 
on large scales, the underlying dark matter density field (for a re- 
view see Meiksin (2007)). At lower redshifts, the situation is likely 
to be more complicated (e.g. iDave et alj | |2003I) ): the non-linear 
evolution of cosmic structures changes the simple picture above al- 
lowing Lya absorbers to populate a variety of environments from 
the l arge scale st ructure to galaxy groups and underdense regions 
( e.g. iLe Brun et al. (19 96): Penton et a l. (20 02); Rosenberg et aT 
(20031): lLanzettaetal.1 fr9 96): Bowe n et al.l |2p02); McL in et al 
(20ol); iGrogin & GeileJ fl998,) : .Cote et alj ( [2005) : .Putman et al 
(20061)). Furthermore, because of the atmospheric absorption 
of UV-photons, the low redshift Lya absorbe rs can be studie d 
only from space based observatories (Wevm ann et al.| ( Il998h : 
iTripD et al.l j2002h ) on a limited number of lines-of-sight mak- 
ing the results potentially affected by cosmic variance and/or small 
number statistics. The cross-correlation function between low red- 
shift galaxies and Lya absorbers is the cleanest statistic for quanti- 
fying the relation between the two populations and has been in- 
vestigated recently bot h observatio n ally a n d using some hydro - 
dyna mical simulation s jChen et alj j2005l) : iRvan-Weberl ( |2006|) : 
.Wilman et alj j2007l) ), with somewhat contradictory findings. 
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RW06 using t he HI Parkes All Sky S urvey (HIPASS) data set 
jMevej i l2003h : IWong, O. I. et alj ilW^ ) has found a puzzling re- 
sult: the galaxy-absorber cross-correlation signal is stronger than 
the galaxy auto-correlation on scales 1-10 Mpc. Earlier stud- 
ies, based, however, on a li mited sample of 16 Lyg lines-of-sight, 
showed the opposite trend ( iMorris & Jannuzil ( |2006|) ). The RW06 
result is not well reproduced either observationally or theoreti- 
cally by Wilman et al who relied on a different data set 

orris & Jan nuzi ( 2006)) an d considere d a sin gle hydrodynami- 
Chen et al. bOOSh se em to be more 
' ' 20071) . However, it is 



Wilman et al. 



cal simulation. The results of 
consistent with the findings of 
worth stressing that while the RW06 galaxy sample includes low 
redshift objects the other two have been obtained from magnitude 
limited catalogs at higher redshifts. 

In this paper we compute the auto and cross-correlation func- 
tions of more than 6000 Lya absorbers over ~ 1000 independent 
lines-of-sight and ~ 5000 mock galaxies extracted from the z — 
outputs of three different high-resolution hydrodynamical simula- 
tions of a ACDM universe in order to better investigate the issues 
above. 

In Section |2] we present the numerical experiments and de- 
scribe the samples of simulated galaxies and Lya absorbers. The 
details of the auto and cross-correlation analyses are described in 
Section[3] The correlation analysis of the mock samples of galaxies 
and absorbers is performed in real-space (Section |4ll and redshift 
space (Section |5}. The results are then summarized and discussed 
in Sections |6]and [7] 



2 HYDRODYNAMICAL SIMULATIONS AND MOCK 
SAMPLES 

We use a set of three hydrodynamical simulations run with 
GADGET-2 and its new fastest version GADGET-3, a parallel tree 
Smoothed Particle Hydrodynamics (SPH) co de that is based on the 
conservative 'entropy -formulation' of SPH jSpringel & Hemguistl 
l2002l : ISpringeill2005i) . The simulations cover a cosmological vol- 
ume (with periodic boundary conditions) filled with an equal num- 
ber of dark matter and gas particles. Radiative cooling and heating 
processes are followed for a primor dial mix of hydrogen and he- 
lium f ollowing the implementation of lKatz. Weinberg & Hemguisll 
( Il996l) . We assume a mean Ultra Violet Background (UVB) pro - 
duced by quasars and galaxies as given bv lHaardt & Madad fl993) . 
with the heating rates multiplied by a factor 3.3 in order to better 
fit observational constraints on the temperature evolution of the In- 
tergalactic Medium (IGM) at high redshift. Multiplying the heating 
rates by this factor (chosen empirically) results in a larger IGM tem- 
perature at the mean density which cannot be reached by the stan- 
dard hydrodynamical code but aims at mimicking, at least in a phe- 
nomenological way, the non-e quilibrium io n izatio n effects around 
reionization (see for example iBolton et alj ( l2007h ). The star for- 
mation criterion for one of the simulations (No Winds - NW) very 
simply converts all gas particles whose temperature falls below 10^ 
K and whose density contrast is larger than 1000 into (coUisionless) 
star particles, while for other two simulations with strong galactic 
winds (Strong Winds ~ SW and Extreme Strong Winds - ESW) a 
multiphase star formation criterion is used. 

The implementation of galactic winds is described in 
ISpringel & HernquistI J2003h but we summarize here the main fea- 
tures. Basically, the wind mass-loss rate Mw is assumed to be pro- 
portional to the star formation rate, and the wind carries a fixed 
fraction x of the supernova (SN) energy. Gas particles are stochas- 



Simulations 


Run 


vw (km/s) 


X 


V 




Iw (kpc) 


NW - - - - 


SW 


484 


1 


2 


0.1 


20 


ESW 


484 


2 


4 


0.025 


60 



Table 1. Main parameters of the simulations. NW (No Winds) uses the 
quick option for the star formation criterion that converts all the gas par- 
ticles below 10^ K and above 5 = 1000 into stars. SW (Strong Winds) and 
ESW (Extremely Strong Winds) models use the default multiphase star for- 
mation criterion. The density = Sw Pth denotes the threshold density 
for the decoupling of the hydrodynamic force, and Iw indicates the wind 
free travel length. 



tically selected and become part of a blowing wind, then they are 
decoupled from the hydrodynamics for a given period of time or till 
they reach a given overdensity threshold (in units of pt^ which is 
the overdensity threshold for star formation) in order to effectively 
travel to less dense regions. Thus, four parameters fully specify the 
wind model: the wind efficiency rj, the wind energy fraction x, the 
wind free travel length 1^ and the wind free travel density factor 5^ . 
The first two parameters determine the wind velocity through 
the following equations: 



and 



Mu, = rjM^, 



1-2 



(1) 



(2) 



from which one can compute the maximum allowed time of the 
decoupling tdec = Iw/vw The parameter 1^ has been introduced 
in order to prevent a gas particle from getting trapped into the po- 
tential well of the virialized halo and in order to effectively escape 
from the Inter Stellar Medium (ISM), reach the low density IGM 
and pollute it with metals. We used si milar values to those tha t have 
been adopted by recent studies (e.g. iNagamiiie et. al.l ( 12007 1) that 
found that the outcome of the simulation is relatively insensitive 
to the choice of this parameter. We note that this wind implemen- 
tation is different from the momentum-driven implementation of 
lOppenheimer & Davd ( 12006 ). which seems to better fit statistics of 
CIV absorption in the high-redshift universe. 

Throughout h indicates the Hubble constant at the present 
epoch, .ffo in units of 100 km s^^ Mpc^^. The cosmological 
model corresponds to a 'fiducial' ACDM Universe with Qra = 
0.26, Qa = 0.74, = 0.0463, n, = 0. 95, Hp = 72 k m 
s"^ Mpc"^ and ag = 0.85 (the B2 series of I Viel et al.l ( l2004h ). 
These parameters provide a good fit to the statistical properties of 
transmitted Lya flux at 2 > 2. We use 2 x 400^ dark matter and 
gas particles in a volume of size 60 h^^ Mpc box and the simula- 
tions are evolved down to 2: = 0. The gravitational softening is set 
to 5 h"^ kpc in comoving units for all the particles. The mass per 
gas pa rticle is about 4.3 x 1 0^ Mq which is a factor ~ 5 better than 
that of IWilman et"ai] ilOOj) . 

These three simulations offer us the opportunity to investigate 
the galaxy-IGM interplay at z — taking into account the role of 
different amount of feedback in the form of galactic winds and the 
role of two different criteria of star formation. Note that similar in- 
vestigations using the same hydrodynamical code and focussing on 
the properties of neutral hydrogen around Damped Lya systems 
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have been performed by iNagamine et. ah In Figure [T] we 

present a qualitative view of the neutral hydrogen overdensity in a 
slice of thickness 6 comovine Mpc ft"^ for the ESW run. We note 
a clear tendency for neutral hydrogen to avoid hot environments, 
where the neutral fraction is lower. The HI distribution in the NW 
and SW simulations it is almost identical on the scale of the plot, 
and therefore are not shown here. Differences can only be spotted 
on scales smaller than 0.5 comoving Mpc h~^in which compact 
knots of neutral hydrogen are seen in the ESW that are not present 
in the NW simulation, since the simplified star formation criterion 
of this latter converts cold gas into coUisionless stars. We will ad- 
dress the differences between the simulations in a quantitative way 
in the following sections. 



Mock Galaxy Samples 





"gal 


Mm- 




\ivii3fvr; 




Wind 


Gnw 


4980 


8.0 


3160 


24.6 


4.7 


NW 


HG 


2480 


19 


3160 


53.4 


10.9 


NW 


LG 


2500 


8.0 


19 


11.4 


2.0 


NW 


Gsw 


4980 


8.6 


3128 


25.6 


2.1 


SW 


Gesw 


4980 


8.6 


3100 


25.4 


1.9 


ESW 



Table 2. Mock Galaxy Samples. Column 1 : Sample name. Column 2: Num- 
ber of mock galaxies. Column 3. Minimum dai'k halo mass. Column 4. Max- 
imum dark halo mass. Column 5. Geometric mean dark halo mass. Column 
6. Geometric mean baryonic mass. Column 7. Wind Model. All masses are 
in 10^°Mq ft-l units. 



2.1 Mock galaxies 

In the simulation we assume a one-to-one correspondence between 
gas-rich galaxies and their dark matter halo hosts. We extract ha- 
los using a friend-of-friend algorithm with a linking length which 
is 0.2 times the mean interparticle separation and consider only 
identified haloes in the mass range [8 x 10^", lO^^ '^ M© /i"^]. 
The lower limit is set (conservatively) by the numerical resolution 
while the upper limit avoids including large halos associated with 
groups and clusters, rather than single galaxies. However, we have 
checked that including the few halos larger than 10^^'^ M© does 
not affect the results presented in this work. The geometric mean 
mass of the haloes is ~ 2.46 x 10^^ M0 h^^, to be compared with 
a mean mass 10^^ M© associated to dar k matter halos host- 
ing HIPASS galaxies (RW06. IMo et al.1 ( |2005|) ). The space density 
of these mock galaxies (0.0023 per cubic Mpc comoving) is 
si milar to that of HI PASS galaxies in the volume limited sample 
of iMever et ai] bOOV ) [M07] (~ 0.003 per cubic Mpc h'^). This 
sample contains all galaxies within 30 Mpc h^^ and HI mass above 
10'' °'^ M0 corr esponding t o a ha lo mass of ^ lO" Mq 
as inferred from the IMo et al. I i2005h model, i.e. similar to our 
lower mass cut off. As we shall see in Section |4] the spatial two- 
point correlation function of these mock galaxies matches that of 
the HIPASS objects, hence fulfilling the main requirement of our 
analysis. 

The mock galaxies extracted from the three simulations are 
hosted in the same dark matter haloes that, however, have a differ- 
ent baryon (gas-l-star) content. The baryon mass in the mock galax- 
ies is affected by galactic winds and star formation processes. The 
mean baryonic mass measured in the NW, SW and ESW simula- 
tions is respectively 4.7, 2.1 and 1.9 x 10^" Mq , thus in- 
dicating that galactic winds are quite effective in blowing baryons 
out of dark halos. The star formation mechanism also plays a role: 
the mean stellar mass of 3.0 X 10^° M© h-^ in the NW simula- 
tion decreases to 0.6 and 0.5 x 10^° Mq h'^ in the SW and ESW 
experiments that adopt the multiphase criterion. 

To better investigate the dependence of the spatial correlation 
on the galaxy mass we have divided, for the NW case only, the 
mock galaxy sample by mass in two subsets. The characteristics of 
all mock galaxy samples considered in this paper are summarized 
in Table 2. 

Finally, to compute the correlation properties of the mock 
galaxies we have generated a random galaxy sample by randomly 
positioning 5 x 10^ objects in the simulation volume. 



2.2 Mock Lya absorbers 

The computational box was pierced with 999 straight lines running 
parallel to the three Cartesian axes. Three sets of 333 mock Lya 
absorption spectra along each axis were simulated and analyzed, 
both in real and redshift-space, to measure the position of each 
Lya line and the column density of the associated HI absorber. 
In this work we only consider weak Lya absorbers with column 
densities in the range 12.41 s; log(7VHi/cm"^) < 14.81 to match 
the characteristics of the RW06 sample. 

The total number of absorbers increases slightly in presence 
of winds, while their average column density decreases, as shown 
in Table 3. However, the differences are small, especially between 
the SW and ESW experiments. The density of Lya absorbers along 
the line-of-sight in the NW simulation 10^'^ km~^s) is larger 
than in the RW06 sample (~ 4 x 10~* km^^s). 

To investigate the significance of this mismatch we have com- 
puted the number of Lya absorbers in our mock spectra, per unit 
redshift and column density in e ach of the three sim ulations and 
compared it with that measured bv lPenton et al.l ( |2004|) in the Space 
Telescope Imaging Spectrograph (STIS) QSO spectra. The results 
are shown in fig. [2] The solid, red curve refers to the NW simula- 
tion. The short-dashed green and the dot-dashed blue curves repre- 
sent the Lya lines in the SW and ESW runs, respectively. The dis- 
tribution of the absorbers is robust to th e presence of galacti c winds. 
When compared to the STIS data of IPenton et alj j2004h (long- 
dashed black curve) we note that the number of absorbers predicted 
by the simulation is larger than the observed ones over most of the 
TVhi range sampled by RW06 (indicated by the two vertical dotted 
lines). The difference between models and data, however, is well 
within observational errors of~ ldexforlog(iVHi/cm"^) ^ 14.5 
jPenton et al] j2004h ). Since we expect that similar observational 
errors for RW06 absorbers, we conclude that there is no significant 
difference in the number density of mock and RW06 Lya lines. 

To investigate the dependence of the clustering properties on 
the absorber column density we have set a column density threshold 
TVhi ~ 10^^'^'* cm'^^which divides the sample in two equally large 
subsets and sorted all mock absorbers in the NW simulation by 
column density. 

The main characteristics of each mock absorber sample are 
listed in Table 3. Moreover, since these mock sample contains many 
more spectra than in the real case, we have also extracted several 
absorbers' sub-samples of 27 lines-of-sights to mimic the RW06 
sample and assess the sampling noise. 
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Figure 1. The spatial distribution of the neuti'al hydrogen overdensity in a slice of thickness 6 Mpc (comoving) around the largest cluster in the simulation 
box (bottom-right part of the panel) extracted from the ESW simulation. White dots are drawn at the position of dark matter halos. Their size is proportional 
to the halo mass. 



Mock Absorber Samples 



Sample 


NAbs 


A^HIMin 


A'"HIMax 


Wind 


A NW 


6239 


12.41 


14.81 


NW 


HA 


1917 


13.24 


14.81 


NW 


LA 


4322 


12.41 


13.24 


NW 


Asw 


6444 


12.41 


14.81 


SW 


Aesw 


6445 


12.41 


14.81 


ESW 



Table 3. Absorber Samples. Column 1: Sample name. Column 2: Number 
of mock LyQ absorbers. Column 3. Minimum column density. Column 4. 
Maximum column density. Column 5. Wind model. All column densities 
are in log{cm^'^) units. 



Finally, to compute the two-point spatial correlation func- 
tions, we have generated random absorber samples by randomly 
positioning 50 Lya absorption lines along the same 999 lines-of- 
sight used for the mock Lya absorption spectra. We verified that 
the estimation of the correlation function does not change signifi- 
cantly if, instead, we consider 999 randomly chosen lines-of-sight 
for the random absorber samples. We note that 50 lines per spec- 
tra represents a good compromise between accuracy and comput- 
ing time since doubling the number of random absorbers does not 
modify our estimates of ^. 



-10 



12 




HI Ly-a lines [NW] ^ 
HI Ly-a lines [SW] 
HI Ly-a lines [ESW] 
Penton Stocke & Shull 2004 



1 4 

Log N„, (cm-2) 



Figure 2. Number of Lyo absorbers per unit redshift and column density. 
Sohd red hne: NW simulation. Dashed green lin e: SW simulation. Dot- 
dashed blue line: ESW simulation. Dashed black: jPenton et alj2004h best 
fit. 
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3 CORRELATION ESTIMATORS 

In this work we use the lOavis & Peebles! ( Il983l) estimator to com- 
pute the galaxy-absorber cross-correlation function both in real 
and redshift-space, £,{Tp, tt) as: 

AG(rp,7r) URG 



(3) 



RG{rp,n) UAG 

where AG{rp,n) is the number of mock absorber-galaxy pairs 
with projected separation, rp in the range [vp — 5rp/2,rp + 
Srp/2], and separation along the line-of-sight, tt, in the range 
[tt — 5TT/2,n + 5tt/2]. RG{rp, n) is the number of pairs consisting 
of a random absorber and a mock galaxy. In both axes the binning 
Sa and Sn is set at 0.39 /i^qo Mpc, i.e. 4 times wider than in RW06. 
The pair counts are divided by the total number of random-galaxy 
pairs urg, and galaxy-absorber pairs, uag- The separations rp 
and TT between two objects are computed from th eir recession ve- 
locities V,; and Vj according to dFisheret all 1994) : 

l-s s.s , 



^^o|l| 



s • s 

^0 



where 1 = (vi +\2)/2 and s = Vi — V2. The estimator ^ is evalu- 
ated in the range of separations [0, 50] Mpc /i^^both along rp and 
TT directions. To estimate the galaxy-absorber correlation function 
in redshift-space we have used the distant observer approximation, 
i.e. we have counted the galaxy-absorber pairs in each of the three 
subsets of mock spectra parallel to one Cartesian axis and consid- 
ered only the corresponding component of the peculiar velocity to 
compute the redshift. The rationale behind this choice is to detect 
and average out possible geometrical distortions arising, for exam- 
ple, when lines-of-sights are oriented along Hl-rich gas filaments 
or when a large fraction of mock galaxies belong to some promi- 
nent, anisotropic cosmic structure. 

The galaxy-galaxy and absorber-absorber auto-correlation 
functions are calculated in a similar way, i.e. by counting galaxy- 
galaxy and absorber-absorber rather than galaxy-absorbers pairs. 
The spherical average of ^(a, n) gives the spatial correlation func- 
tion ^(s) where s — \Jr'p + tt^. We also estimate the analogous 
quantity in real-space, ^(r), where r represents the genuine pair 
separation that coincides wit their redshift difference in absence 
of peculiar velocities. In order to compare our result with those of 
RW06 we compute two more quantities. The first one is the pro- 
jected correlation function, H(rp): 



E{rp) = 2 / 
Jo 



C(rp,7r)d7r, 



(5) 



where TVmax = 50 Mpc h~ . 

The second one is the absorber auto-correlation along indi 
vidual lines-of-sight, 

AA{-!V) UAR , 



(6) 



AR{n) riAA 

where AA{-k) is the number of mock absorber pairs with separa- 
tion TT along the line-of-sight and AR(-k) is the number of random 
absorber pairs. 

The uncertainties in the cross and auto-correlation functions 
of the mock samples are computed using the bootstrap resam- 
pling technique. For large, independent datasets bootstrap errors 
are equivalent to uncertainties calculated using the jackknife re- 
sampling, as in RW06. The uncertainty is computed in each (fp, tt) 
bin as 



2 '^1 = 1 (^^ ) 



iV- 



(7) 



where the subscript i identifies the bin, j refer the sample and is 
the average correlation function computed over the N bootstrapped 
samples. In this work TV = 50 which provide us with a robust error 
estimate (increasing to 350 modifies errors by < 2%). 

This error estimate assumes that the covariance matrix of the 
data is diagonal, i.e. that the values of ^(cr, tt) in different bins are 
not independent, which is known not to be the case. However, our 
simple way of estimating the uncertainties avoids the complication 
of dealing with a large covarian ce matrix, while prov iding an unbi- 
ased estimate of the real errors faawkins et"ai]|2003h . 



4 REAL-SPACE ANALYSIS 

In this analysis we ignore peculiar velocities when we use Eq. |4] 
to estimate rp and tt from redshifts. In Fig. [3] we show the real- 
space auto-correlation function of the mock galaxies in the Gnw 
sample (black dots). Errorbars represent l-cr bootstrap uncertain- 
ties. The autocorrelation of mock galaxies is shown together with 
that of HIPASS galaxies, indicated by the dashed line which rep- 
resents the power-law best fit to the ^(r) in the volume-limited 
sub-sample of galaxies extracted from the HIPASS catalog by MOT. 
This power-law has a slope 7gg = 1.5 ± 1 and correlation length 
ro,gg ~ 3.2 ± 1.4 Mpc . The two functions agree, within the 
errors, below 10 Mpc , since the power-law fit to the correla- 
tion function of our mock galaxies in the range [1,10] Mpc ft^^has 
-igg = 1.46 ± 0.03 and ro,gg = 3.06 ± 0.15 Mpc We have 
considered the M07 result since it is based on a sub-catalog that is 
volume limited, like our mock samples but it is worth noticing that 
the RW06 fit obtained using the full, flux limited HIPASS sample 
is fully consistent with the M07 result and, therefore, with our fit 
too. 

The correlation signal of the mock galaxies suddenly drops at 
separations smaller than 0.4 Mpc . On the contrary, the galaxy 
correlation function of RW06 monotonically increases when re- 
ducing the pair separation. Including the few mock halos larger 
than IO^'^'^Mq sample does not modify significantly this small- 
scale trend. This small-scale mismatch as an artifact deriving from 
the fact that, in the simulation, we do not resolve galaxy-size 
sub-structures within the large cluster-size halos that, if present, 
would significantly contribute to the correlation signal at sub- 
Mpc h,"'^ scales. Indeed, when we run the Friends-of-Friends algo- 
rithm to identify halos using a smaller linking length of 0. 1 times 
the mean inter-particle spacing, the small scale flattening disap- 
pears and the power-law behavior is restored below 0.3 Mpc h^^ . 

Mock absorbers are significantly less self-clustered than 
galaxies: their autocorrelation functi on (blue squares) i s factor of 
~ 10 below that of galaxies (see Dob rzvcki et alj ( l2002h ). We can- 
not compare this result with observational data directly, since the 
observed Lya absorbers are too sparse. However, RW06 was able 
to compute their correlation along each line-of-sight and we com- 
pare this result with the theoretical predictions in the next Section. 

The red triangles show the mock galaxy-absorber cross- 
correlation function of the Anw+Gnw samples which is signif- 
icantly weaker than the galaxy auto-correlation. This result is at 
variance with that of RW06 who find that the cross-correlation 
function of HIPASS galaxies and Lyo absorbers (dot-dashed curve 
in Figure [3]l in the [1, 10] Mpc ft^^ range is best fitted with a 
power-law slope 7ag = 1.9 ± 0.3 and correlation length ro^ag ~ 
7.2 ± 1.4 Mpc h~^, significantly larger than that of the galaxy 
auto-correlation function. When we fit the cross-correlation func- 
tion of the mock data in the same range of separations we find 
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Figure 3. The real-space two-point correlation functions of tlie mock 
galaxy and absorber samples. Black dots: galaxy auto-coirelation func- 
tion in the G^W'Sample. Blue squares: absorber autocorrelation function 
in the ANW-sample. Red triangles: galaxy-absorber autocorrelation func- 
tion in the Anw+Gnw sample. The size of the bars shows l-cr bootstrap 
resampling uncertainties. Black dashed curve: best fit to the galaxy-galaxy 
congelation function in the HIPASS volume-limited sample of M07. Black 
dot-dashed curve: RW06 best fit to the HIPASS galaxy- hya absorbers 
cross-correlation function. 



Jag = 1.29 ± 0.03 and correlation length ro,ag = 1.44 ± 0.08 
Mpc h'^. 

RW06 pointed out that the cross-correlation signal increases 
with the column density of the absorber. We find the same trend 
in the simulation. We show in Fig. |4l the cross-correlation signal 
increases when we restrict our analysis to strong absorbers of the 
HA sample. On the contrary, massive mock galaxies do not seem 
to be significantly more or less correlated to Lya absorbers than 
smaller galaxies. In fact, we find that the cross-correlation signal is 
almost independent of galaxy mass. 

Strong galactic winds can blow gas out of galaxy-size halos 
and therefore could suppress the cross-correlation signal on sub- 
Mpc scales. To quantify the effect we have computed the galaxy- 
absorber correlation functions in the SW and ESW simulations 
and compared them with that of the NW experiment. The results 
are shown in Fig. [5] The red triangles with errorbars represent the 
same cross-correlation function of the Anw+Gnw sample shown 
in Fig. |4] and refer to the case of no winds. The effect of includ- 
ing the effect of strong winds is illustrated by the blue dashed and 
solid black curves that refer to the SW and ESW simulations, re- 
spectively. Even adopting extreme prescriptions for galactic winds, 
the effect on the galaxy-absorber correlation function is very small 
and, as expected, is significant only at separations < 0.3 Mpc 
where fewer galaxy-absorber pairs are found with respect to the 
NW case. This is not surprising, considering the free travel length 
Iw adopted in the models. We find no significant differences be- 
tween the SW and ESW experiments, which illustrates the robust- 
ness of the cross-correlation signal on scales larger than /n, to the 
scheme adopted to simulate galactic winds. 



Figure 4. The absorber-galaxy cross-correlation function in various mock 
subsamples. Red triangles with en'orbars: Anw+Gnw sample. Black soHd 
line: HG-hAnw- Dashed red: LG-i-Anw- Dotted blue: Gnw+HA. Dot- 
dashed green: Gnw+LA. Long-dashed magenta: LG-l-HA. 
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Figure 5. Effect of galactic winds on the cross-correlation function. Red 
triangles: Anw+Gnw samples in the NW simulation. Dashed blue curve: 
SW simulation. Solid, black curve: ESW simulation. The size of the bars 
shows l-cr bootstrap resampling errors. 



5 REDSHIFT-SPACE ANALYSIS 

In Section |4] we have shown that hydrodynamical simulations 
do not reproduce the RW06 result. On the contrary, the galaxy- 
absorber correlation function is significantly weaker than the 
galaxy autocorrelation function. The previous analysis, however, 
has been performed in real-space ignoring peculiar velocities that 
may bias the correlation analysis. Moreover we have considered a 
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number of spectra much larger than that of RW06. Therefore, we 
must account for the possibihty is that the mismatch between hy- 
drodynamical simulations and RW06 is not genuine but derives, in- 
stead, from redshift-space distortions and sparse HI sampling that, 
if not properly accounted for, may affect the cross-correlation anal- 
ysis. In an attempt to account for both types of errors we repeat the 
correlation analysis using more realistic mock catalogs in which 
redshifts are used as distance indicators and only 27 lines-of-sight 
are taken to mimic the RW06 data set. To investigate the two ef- 
fects separately, we first perform a redshift-space analysis of the 
whole Anw+Gnw sample and then we consider sub-samples of 
27 lines-of-sight. 

In Figure [6] the autocorrelation function of the mock galax- 
ies in the Gnw sample, ^gg{rp, tt), is plotted on the {vp, tv) plane. 
Contours are drawn at iso-correlation levels of 2,1,0.5,0.25. The 
distortions along the tt axis induced by small scale incoherent mo- 
tions within virialized structures (the so called fingers-of-god) can 
be seen at separations rp ^ 2Mpc h^^ extending out to tt = 6 
Mpc h~^. A similar distortion pattern is seen in the correlation 
function of HIPASS galaxies (Fig. 2 of RW06). In that case a sec- 
ond, independent, distortion pattern along the Vp axis, is detected 
at separations rp > 4Mpc /i^^. The compression of the isodensity 
contours along tt is the signature of large scale coherent motions 
that increase the apparent number of pairs with large separations. 
This second distortion pattern is not visible in Fig.|6l a fact that we 
ascribe to the lack of large scale power in our simulations. Indeed, 
our simulations do not account for power on scales larger than 60 
Mpc /i~^which could significantly contribute to the amplitude of 
the bulk motions and thus to the compression of the iso-density 
contours. 

Figure |7] shows the redshift-space cross-correlation function 
^ag{rp,n) of mock absorbers and galaxies in the Anw+Gnw 
sample. The signal is significantly weaker than the galaxy auto- 
correlation and the distortion pattern looks very different as no sig- 
nificant elongation is seen along the tt axis. Instead, at large sep- 
arations, the iso-correlation contours are compressed along tt, as 
expected in the presence of coherent motions. The differences be- 
tween ^agivp, tt) and ^gg{rp, tt) reveal that mock Lya absorbers 
and galaxies have different dynamical properties. Galaxies' rela- 
tive velocities are dominated by the incoherent motions, typical of 
virialized structures. Instead, the relative motion of mock Lya ab- 
sorbers and galaxies is more coherent, suggesting that mock ab- 
sorbers are preferentially located in the outskirts of high density 
regions into which they are probably falling. 

Finally, we note that the peak of the cross-correlation function 
is spatially offset from the center. This feature and the general dis- 
tortion pattern of the simulated cross-correlation function are quali- 
tatively similar to that of the cross-correlation function between the 
CHFT galaxies and the HST Quasar Absorption Line Key Project 
Data Release Lya with 13 ^ log(A'Hi/cm~^) < 15 measured by 
W07. On the contrary, the RW06 cross-correlation function is dom- 
inated by a very large finger-of-god distortion. A similar, but less 
prominent, distortion pattern has been seen by Daveetal. (1999) 
and W07 in their numerical experiments. RW06 interpreted this dis- 
tortion as the draining of the gas from low-density regions into col- 
lapsed structure. Although the dynamical interpretation in this case 
is not as simple as in the galaxy-galaxy case, we note that the drain- 
ing mechanism advocated by RW06 would probably lead to coher- 
ent, rather than incoherent motions, which would produce a very 
different distortion pattern. W07 suggested that the finger-of-god 
distortion could be a geometrical effect deriving from observing 
Lya absorbers along lines-of-sights that run along some radially- 




Figure 6. The redshift-space auto-correlation function ^gg (rp , n) of mock 
galaxies in the Gnw sample. The binning is 0.4 Mpc h~^'m both axis. Iso- 
correlation contours are drawn at correlation levels of 2.1,0.5,0.25. Boot- 
strap resampling shows typical pixel variations of the order C()"p,7r) ~ 
0.08 



elongated structure. To check this hypothesis we exploited the dis- 
tant observer approximations and computed the cross-correlation 
function by considering redshift distortions along one Cartesian 
axis at a time. If distortions were purely geometric, i.e. induced 
by a few prominent, anisotropic structures, we would expect to see 
different distortion patterns in the cross-correlation functions com- 
puted along orthogonal axes. If, on the other hand, they were caused 
by random motions within large, spherically symmetric, virialized 
structures like galaxy clusters, we would expect to see fingers-of- 
god type distortions along all axes. Instead, the correlation func- 
tions measured by three, orthogonally-positioned distant observers 
turned out to be very similar and consistent with the one shown 
in Fig. [7] We conclude that neither pure geometrical effects nor 
incoherent motions can alone explain the distortion pattern in the 
^ag{rp, n) of our mock Anw+Gnw samples. 

Small redshift distortions could be amplified by sampling 
Lya absorbers along a limited number of lines-of-sights, as in 
the RW06 case. To quantify the effect of shot noise errors cou- 
pled to dynamical and geometrically-induced distortions, we have 
constructed 30 independent realistic mock Lya sub-samples of 27 
independent lines-of-sights and computed their cross-correlation 
with all mock galaxies of the Gnw sample. In Fig. [8] we show 
^ag{r p,n) computed in four such realistic mock samples. The 
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Figure 7. The redshift-space cross-correlation function ^agirpjir) of 
mock absorbers and galaxies in the Anw+Gnw sample. The binning is 
0.4 Mpc h^^in both axis. Iso-correlation contours are drawn at correlation 
levels of 1,0.5,0.25. Bootstrap resampling shows typical pixel variations of 
the order 5(rp, tt) 0.2 



cross-correlation functions shown in the two upper panels are char- 
acterized by prominent finger-of-god distortions which, in the 
upper-right plot, are similar in amplitude to that measured by 
RW06. This kind of distortion is found in ^ 20% of the mock 
subsamples considered. The fact that we observe fingers-of-god 
distortions along different Cartesian axes suggest that they cannot 
be attributed to the fact that the sample is dominated by a sin- 
gle, prominent, anisotropic structure. Rather, they seem to originate 
from genuine, finger-of~god like, dynamical distortions which be- 
come apparent when a significant fraction of the 27 spectra samples 
some virialized regions. The relevance of sparse sampling variance 
in the cross-correlation analysis is even more evident in the two 
bottom panels of Fig. [8] They show the cross-correlation function 
computed along the same (Z~) axis, as in the top-right panel, but 
use two independent sets of lines-of~sight. Not only the finger- 
of-god distortion disappears but the cross-correlation signal is ei- 
ther very weak (bottom left) or significantly offset from the center 
(bottom-right). 

A more quantitative assessment of sparse sampling errors is 
given in Fig. [9] in which we show the projected absorber-galaxy 
cross correlation function 'E,ag(rp)/rp of the Anw+Gnw sam- 
ple (filled black dots). Small errorbars drawn with solid lines rep- 
resent l-(7 bootstrap resampling errors computed using all 999 



mock absorbers in the A catalog. Large errorbars plotted with 
dashed lines represent the scatter around the mean of the projected 
cross-correlation function computed using the 30 realistic mock 
absorbers' samples consisting of 27 lines-of-sight. The sampling 
noise clearly dominates the error budget and the total error signifi- 
cantly exceeds that of RW06. Filled red squares show the projected 
galaxy-galaxy correlation function with the 1-a bootstrap errors. 
In order to assess the goodness of our error estimate we have com- 
pared the scatter among the 30 catalogs with the bootstrap errors 
computed from N=50 samples. The two errors agree well in the 
range (1,10) Mpc/h, in which boostrap errors are ~ 15 % smaller 
than those shown in Fig. |9] On smaller scales the bootstrap re- 
sampling technique overestimates the errors by factor of ~ 2. The 
autocorrelation signal is higher than the cross-correlation one, con- 
sistently with the real-space analysis. However, the difference is of 
the order of the errors, i.e. the mismatch is about 1-cr at separations 
Tp > 1 Mpc h^^,in the range in which RW06 find that the cross- 
correlation signal is larger than the autocorrelation one. Filled tri- 
angles show the projected autocorrelation function of all absorbers 
in the Anw sample. As anticipated by the real-space analysis, ab- 
sorbers correlate with themselves very weakly. When one accounts 
for sparse sampling their autocorrelation signal is consistent with 
zero. 

RW06 were able to detect the auto-correlation signal of the ab- 
sorbers by measuring their auto-correlation function ^(tt) of eq. |6] 
along individual lines-of-sight. We have repeated that analysis us- 
ing all absorbers in the A sample. The resulting auto-correlation 
function replicates the RW06 result to within l-a. 

Finally, to test the robustness of our results we have computed 
the cross-correlation function, (,ag{rp, n) using the HG, LG, and 
LA sub-samples as well as the mock catalogs extracted from the 
SW and ESW runs. There are no cases in which wer are able to ob- 
tain a galaxy-galaxy autocorrelation signal weaker than the cross- 
correlation one and to reproduce the large finger-of-god distortion 
feature observed by RW06. 



6 CONCLUSIONS 

In this work we have studied the relative spatial distribu- 
tion of galaxies and weak Lya absorbers with 12.41 ^ 
log(AfHi/cm~'^) ^ 14.81 in hydrodynamical simulations and 
compared our results with the analyses of real datasets performed 
by W07 and, mainly, with RW06. Our main conclusions are: 

• The galaxy-absorber two-point cross-correlation function in 
the hydrodynamical simulation is weaker than the galaxy auto- 
correlation function. This result is at variance with that of RW06 
but in qualitative agreement with the analysis of W07. 

• No flattening at small separation is observed in the cross cor- 
relation function of all mock absorbers, unlike in RW06. A small- 
scale flattening is observed, however, when the cross correlation 
analysis is restricted to low density absorbers. 

• The cross correlation signal increases with the column density 
of the absorbers, in agreement with RW06. We find no significant 
dependence on galaxy mass. 

• Galactic winds have a small effect on the absorber and galax- 
ies correlation properties in these models. Using the most ex- 
treme prescription to simulate these winds suppresses the cross- 
correlation signal only at separations < 0.3 Mpc h^^. 

• Absorbers correlate with themselves more weakly than with 
galaxies. Their auto-correlation signal is very weak and consistent 
with that measured by RW06. 
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Figure 8. The redshift-space absorber-galaxy cross-coiTelation function ^ag{rp,iT) for four independent subset of absorbers along 27 lines-of-sight. On 
top each panel is indicated the direction along which absorptions spectra were drawn and the total number of mock absorbers in each sample N^bs ■ All 
cross-correlation functions are computed using the same 4980 mock galaxies in the G sample. 



• Redshift-space distortions alone cannot explain two aspects of 
the differences with the RW06 results. The cross-correlation signal 
is weaker than the galaxy auto-correlation signal. The two-point 
cross-correlation function, (,ag{rp,Tv) does not show a prominent 
finger-of-god type of distortion. The latter looks very prominent 
in the RW06 cross-correlation function but is not seen in the W07 
one. 

• The origin of the finger-of-god distortion cannot be purely 
geometric, i.e induced by the presence of a prominent, anisotropic 
structure in the sample. In this case distant observers taking spectra 
along orthogonal directions would detect different distortion pat- 
terns. We do not see such effect. 

• Fingers-of-god distortions may appear when sampling the in- 
tergalactic gas using a limited number of UV spectra, as in the 
RW06 sample. In this case, they represent genuine dynamical dis- 



tortions that become apparent when a few spectra, that however 
represent a significant fraction of the total, pierce some virialized 
regions. 

• The sampling noise is large. Once accounted for, the differ- 
ence between the simulated galaxy-galaxy and galaxy-absorber 
correlation functions is significant at the ~ 1-cr level only. 



7 DISCUSSION AND PERSPECTIVES 

Modeling the gas distribution in the low redshift universe is a diffi- 
cult task. Numerical experiments use a number of simplifying hy- 
pothesis and approximations that potentially affect our results. The 
main uncertainties are related to the ill-known mechanisms of stel- 
lar feedback and galactic winds for which we have adopted simplis- 
tic phenomenological prescriptions. It is therefore very reassuring 
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Figure 9. Filled dots and solid line: projected absorber-galaxy cross cor- 
relation function 'Eag{rp)/rp. Small errorbars are 1-cr bootstrap resam- 
pling errors. Large errorbars account for sparse sampling variance. Filled 
squares and dashed line: galaxy-galaxy projected autocorrelation function 
^ag{Tp) /rp. Filled triangles and dot-dashed line: absorber-absorber pro- 
jected autocorrelation function "Eaaifp) /^p 



that our results are robust to the star formation criterion and galac- 
tic wind prescriptions adopted. However, since robustness does not 
exclude systematic errors one needs to be aware that the various 
approximations adopted in our numerical model to predict the HI 
distribution at z = may bias our results. 

Our model does not include halos larger that 10^^'^ Mq 
and ignores substructures within virialized halos. While we have 
checked that including large halos does not change our results, 
ignoring galaxy-size halos within groups or clusters may affect 
the outcome of the correlation analysis. Galaxies in strongly clus- 
tered environments significantly contribute to both the auto- and the 
cross-correlation function at small separations. Ignoring their pres- 
ence would artificially decrease the correlation signal, producing a 
flattening in the correlation functions at small separations. We do 
see a flattening but only in the galaxy-autocorrelation function and 
on scales smaller than 0.4 Mpc . The cross-correlation func- 
tion, instead, increases at small separations unlike the one of RW06 
that flattens and we do not reproduce the flattening at separations 
smaller than 1 Mpc h^^ . A flattening of the galaxy-absorber cross- 
correlation fu nction at small sc ales was also seen in the numerical 
simulations of lOave et alj l ll999il that, however, have a limited res- 
olution compared to ours. The fact that we find no flattening in 
the cross-correlation function has two implications. First, ignoring 
sub-clustering within large halos has little impact on our results. 
Second, it seems that there is no characteristic scale for the cosmic 
structures in which Lya absorbers are embedded. 

The RW06 analysis convincingly rules out minihaloes for the 
confinement of weak Lya absorbers. Based on the measured cross- 
correlation strength, RW06 suggest that they are embedded in much 
larger halos with the typical mass of a galaxy-group. This would 
imply a self-clustering of the absorbers comparable or even larger 
than that of galaxies. The fact that, on the contrary, the measured 
absorber self-clustering along the line~of-sight is weak is not re- 



garded by RW06 as a conclusive evidence since redshift distortions 
may artificially dilute the correlation signal. Our numerical exper- 
iments provide a direct estimate for the self-clustering of the ab- 
sorbers which is free of redshift-distortions. The real-space analy- 
sis we have performed indicates that the auto-correlation function 
of the mock absorbers is significantly weaker than that of mock 
galaxies and that, in redshift-space, their self-clustering is con- 
sistent with the RW06 estimates. The outcome of our numerical 
model suggest therefore that in a ACDM universe weak Lya ab- 
sorbers are not embedded in group-size halos. In fact, the asso- 
ciation of weak Lya absorbers with virialized halos is probably 
too naive. The absence of a strong finger-of-god distortions in the 
simulated £,acj{rp,-n) suggest that the neutral hydrogen responsi- 
ble for weak Lya absorption lines is not part of virialized struc- 
tures. Rather, it is probably located in their outskirts, in-falling to- 
wards their central regions. Interestingly, we see a flattening in the 
absorber auto-correlation function at separations < IMpc a 
feature which is also typical of the the w arm-hot intergalactic gas 
according to both nu merical ( lDaveetal.ll2001i) and semi-analytic 
dValageas et"ai]|2002l) predictions. 

Finally, we turn to what we regard as the main result of this 
work. RW06 find that the galaxy-absorber cross-correlation signal 
is significantly larger than the galaxy-galaxy correlation. Our nu- 
merical analysis is not able to reproduce the observation as we 
find that the opposite is true. However, when shot noise errors 
are accounted for, the discrepancy between the auto- and cross- 
correlation signals is of the order of l-tj only. Can we reconcile 
the two results ? Our numerical experiments were performed on 
a rather small box of 60 Mpc which cannot be regarded as 
a fair sample of the universe. In other words our cosmic variance 
is not negligible and should be accounted for in our error budget. 
This would require running numerical simulations in a larger box 
while keeping the same resolution or running several identical sim- 
ulations of different random realizations of the universe. In either 
case the likely outcome would be that of increasing the size of the 
errorbars in fig.|9]and the conclusion would be that, probing the HI 
distribution with 27 lines-of-sight is not sufficient, in a ACDM, 
universe to demonstrate a difference between the self and cross 
clustering of galaxy and Lya absorbers at the level measured by 
RW06. 

The fact that the errorbars in the projected cross-correlation 
function of RW06 are smaller than ours seem to indicate that their 
error estimates are biased low. In section |5] we have shown that 
the bootstrap technique underestimates errors by ~ 15 %, on av- 
erage, at separations ^ IMpc h^^ when the sampling is as sparse 
as in the RW06 case. This bias reflects the fact that absorbers are 
not guaranteed to be independent. It is plausible that this effect is 
even more severe in the RW06 sample since nearly 30% of the ab- 
sorption spectra considered were drawn in the vicinity of the Virgo 
cluster region. We would also expect that these spectra could artifi- 
cially amplify the cross-correlation signal since the Virgo cluster is 
an Hl-rich region. However, surprisingly enough, the excess cross- 
correlation signal is still present when galaxies and absorbers from 
this region are excluded from the analysis (Ryan- Weber, private 
communication) . 

The only way out at this apparent paradox is that the relative 
distribution of galaxies and Lya absorbers in the RW06 sample 
is different from that of the typical cosmic environment, since the 
cross-correlation signal and its variance are significantly different 
from their average values. This despite the fact that in our cosmic 
neighborhood the most prominent structures are anisotropically lo- 
cated along the Super-Galactic plane, rather than being homoge- 
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neously distributed. We see two possible ways to check the validity 
of this hypothesis. One is to resort to the so called constrained hy- 
drodynamical experiments designed to match the actual gas distri- 
bution in our local universe ( Kravtsov, Klypin &. Ho ffman 20021: 
Klvpin. H offman. K ravtzov. & Gottlobej boOsT lYoshikawa et al.l 
2005l : IVieletal.,,206ir Currently available simulations, however, 
are of little use as their constraints are either too weak, as they refer 
to scales larger than 5 Mpc (Gaussian), or too local, as they are 
effective out to distances of ~ 15 Mpc h^^ , i.e. within our local 
Supercluster. The second possibility, which looks more promising, 
is to improve the sampling of the HI distribution either through 
Lya absorption lines in the UV absorption spectra or through the 
X-ray lines of highly ionized metals, like OVII. The latter is ex- 
pected to trace the Warm Hot Intergalactic Medium (WHIM) in 
density-temperature environment similar to that in which the weak 
Lyof absorbers can b e found. With t his respect, proposed X-ray 
satellites like EDGE jPiro et alj[2003) are particularly interesting, 
as they could observe the WHIM in emission, which would allow 
one to trace the three dimensional gas distribution rather than prob- 
ing it in ID along a few lines-of-sight. 

It is worth stressing that the present tension between model 
and data could be a signature of the fact that hydrodynamical sim- 
ulations are still missing physical inputs able to reproduce the ob- 
servations. However, if the mismatch between theory and obser- 
vations is confirmed, which probably requires both better obser- 
vational data and better control over systematics in the numerical 
models, the RW06 results could constitute an interesting challenge 
to the ACDM paradigm, si milar, and perh aps related to the absence 
of dwarf galaxies in voids jPeebleslbOOTh . 
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